Supercluster Magnetic Fields and Anisotropy of Cosmic Rays above 10 19 eV 
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We predict energy spectra and angular distributions of nucleons above 10 19 eV that originate from 
sources distributed in the Local Supercluster, which is also supposed to contain a large scale magnetic 
field of strength ~ 0.05 — 0.5 /j,G. We show that this model can explain all present-day features of 
ultra-high energy cosmic rays, at least for field strengths close to 0.5 fiG. The large-scale anisotropy 
and the clustering predicted by this scenario will allow strong discrimination against other models 
with next generation experiments. 
PACS numbers: 98.70.Sa, 98.62.En 



Ryu & Biermann Q have recently argued that the ob- 
servational upper limit on the strength B ims of an extra- 
galactic magnetic field (EGMF), obtained from Fara- 
day rotation observations of distant sources Ppl], reads: 
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1 fj,G, for fields contained inside the cosmolog- 



ical large-scale structure, such as the Local Superclus- 
ter. Such a strong magnetic field would have profound 
consequences on the propagation of charged ultra-high 
energy cosmic rays (UHECRs) with energy E > 10 EeV 
(1 EeV= 10 18 eV) §-§. In particular, for B Ims ~ 0.1 fiG, 
charged UHECRs with energies up to ~ 100 EeV would 
diffuse, while UHECRs of higher energies would propa- 
gate in nearly straight lines. This would allow to repro- 
duce very nicely the observed energy spectrum of UHE- 
CRs for a single injection spectrum oc E~ 2A [HQ. More- 
over, the associated angular deflection might explain why 
no astrophysical counterpart within ~ 50 Mpc could be 
associated to the highest energy events. The existence 
of a magnetic field of strength ~ 0.1/iG could thus rec- 
onciliate models in which the UHECRs are protons ac- 
celerated in conventional astrophysical sources with 
present-day observations. 

The most recent results of the AGASA experiment || 
provide tight constraints on these ideas and more gener- 
ally on any scenario of UHECR origin. In particular, this 
experiment now reports the detection of 7 events above 
100 EeV, scattered across half the sky, with no obvious 
association with the Supergalactic plane. A naive one- 
source model, as considered in Refs. which predicts 
an angular size of the image ~ 15° for E > 100 EeV, 
is thus excluded. However, whereas AGASA found no 
significant large scale anisotropy, the data indicate sig- 
nificant small scale clustering. In this letter we demon- 
strate that both the observed spectrum and the angular 



distribution can be explained by a diffuse distribution of 
sources with a density proportional to the matter den- 
sity in the Local Supercluster, provided this structure is 
permeated by magnetic fields of strength B lmB ~ 0.5 /zG, 
with power concentrated on ~ Mpc scales. In this sce- 
nario, the large scale isotropy observed by AGASA is 
explained by diffusion, whereas the small scale clustering 
is due to magnetic focusing in the magnetic field struc- 
ture. We discuss future observational tests of this sce- 
nario against other models, especially with regard to the 
strong increase in UHECR statistics anticipated from the 
Pierre Auger Observatory (To) . 

Energy spectrum and angular images. During their 
propagation, UHECR nucleons lose energy by pion pro- 
duction (for E > 50 EeV) , and pair production (protons 
only), on the cosmic microwave background. Charged 
UHECRs also acquire stochastic deflection, and hence 
time delays with respect to straight line propagation, 
in random magnetic fields. Detailed predictions for the 
energy spectrum and angular distribution of UHECRs 
propagating in a magnetic field can only be made through 
numerical Monte-Carlo simulations p Jlljjr^j , in order to 
take into account stochastic energy losses, stochastic de- 
flection, and effects of an anisotropic geometry. In order 
to correctly reproduce the small scale angular distribu- 
tion of events in a given realization of the magnetic field, 
one needs to accomodate a very small solid angle, that 
represents the detector as seen from the source, with a 
reasonable consumption of CPU time. Each of the simu- 
lations presented below typically requires one to several 
weeks of CPU time on DEC ALPHA 500 computers. De- 
tails on our method can be found in Refs. IJUQ. 

We assume that the sources of UHECRs are distributed 
according to the matter density in the Local Superclus- 
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ter, following a pancake profile with scale height of 5 Mpc 
and scale length of 20 Mpc; the observer is located 20 Mpc 
away from the center of the density profile (the Virgo 
cluster), and within 2 Mpc from the middle plane. Fur- 
thermore, we assume that no sources are present within 
2 Mpc from the observer. This represents a reason- 
able modeling of our location in the Local Supercluster 
and of its shape. We also assume that the Local Su- 
percluster is permeated by a random magnetic field of 
strength 0.05 — 0.5 /j>G, whose power spectrum follows a 
Kolmogorov law {B{k) 2 ) oc k nB , with n B = -11/3. The 
largest eddy, defined as the scale over which the phase 
of the magnetic field changes by 2ir, is 10 Mpc, and the 
corresponding smallest eddy is 1 Mpc. This latter value 
is limited by resolution, but it does not influence our re- 
sults, as long as the power is concentrated on the largest 
scales. Some of the dependencies on the strength of the 
magnetic field, the power law index ns of the magnetic 
power spectrum, the small scale cut-off, and the offset 
of the observer position from the middle-plane, are dis- 
cussed below. 
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FIG. 1. Best fit to the combined data above 10 EeV from 
the Haverah Park Jl|, Fly's Eye Q, and AGASA § ex- 
periments (error bars) of the spectra predicted by the diffuse 
source model explained in the text. The thick histogram is 
for -Brms = 0.5 j-i G, with the observer 2 Mpc above the Super- 
galactic plane, and the thin histogram is for Brms = 0.05 /i G, 
with the observer in the plane center. In both cases the spec- 
tra were averaged over 4 magnetic field realizations with 20000 
particles each. 

For a diffuse source distribution, any value of the mag- 
netic field strength, between ~ 0.05 fj,G and ~ 0.5 fj,G 
can provide a reasonable fit to the observed energy spec- 
trum, for a single power-law injection spectrum oc E~ 2A . 
This is illustrated in Fig. |l|, for the above two extreme 
values of B rms . We note that the observed energy spec- 
trum is well approximated by a power-law oc E~ 2,7 in 
the range E < 100 EeV ||. The difference in indices be- 
tween the injected and the propagated spectrum in Fig. [l] 
results from the diffusion of UHECRs combined with en- 



ergy losses 0-0] ; in our present case, UHECRs with en- 
ergies greater than ~ 100 — 1000 EeV propagate in nearly 
straight lines, in which case the energy spectrum is un- 
affected by magnetic deflection, while UHECRs of lower 
energies diffuse. The interplay between propagation in 
the magnetic field and energy losses leads to a spectrum 
that is softer than the injection spectrum in the diffu- 
sive regime, but not in the rectilinear regime. This leads 
to reasonable fits of the UHECR spectrum above ~ 10 
EeV. We emphasize that limiting the injection to an up- 
per cut-off ~ 1000 EeV would not change the spectral 
shape below ~ 500 EeV. 

60. < E/EeV < 9000. 
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FIG. 2. The angular distribution in Galactic coordinates 
of events above 60 EeV, averaged over 4 magnetic field re- 
alizations with 20000 particles each for the scenario with 
-Brms = 0.05 n G, corresponding to the thin histogram in 
Fig. |l|. The grey scale represents the integral flux per solid 
angle. The solid line marks the Supergalactic plane. 
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FIG. 3. Same as Fig. |2|, but for the scenario with 
Brms = 0.5 n G, corresponding to the thick histogram in 
Fig. | 

The angular distributions associated with these simu- 
lations are shown in Figs. ^, and ||, plotted in Galactic 
coordinates, for E > 60 EeV. These images are averaged 
over different spatial realizations of the magnetic field in- 
side the Local Supercluster. For the case of strong fields 
and large coherence lengths, images corresponding to dif- 
ferent realizations are very different from each other due 
to cosmic variance, but cover more than half of the sky, 
consistent with the isotropy observed by AGASA. 

For B rms ~ 0.05 ^G, shown in Fig. |, UHECR arrival 
directions are strongly clustered around the center of the 
Local Supercluster, i.e. approximately around the Virgo 
cluster, at I ~ 282°, b ~ +75°. Since AGASA recorded 
no event out of 47 for E > 40 EeV within ~ 15° of 
this point, which is furthermore located near the peak of 
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the AGASA exposure curve, one would exclude this sce- 
nario to a high degree of confidence. However, we note 
that radio-galaxies in the Supercluster seem to distribute 
rather uniformly with radius, i.e. there is no preferred 
center [H . If the UHECR sources were to follow radio- 
galaxies, the sky distribution would be more isotropic, 
and further statistics would be needed to discriminate 
this scenario (see also below). 

For a field strength B lms ~ 0.5 fj,G, the correlation with 
the Supergalactic plane disappears, as shown in Fig. [|. 
This effect is not trivial, as all UHECR sources in our sce- 
nario are located within the Supergalactic plane. Note 
as well that the resulting angular image in Fig. |3| reveals 
spikes of UHECR events even though the source is dif- 
fuse. These spikes are produced by magnetic lensing and 
give rise to small scale clustering. As we quantify be- 
low, we find that the amplitude of clustering depends on 
the spectrum of magnetic inhomogeneities; namely, for a 
scale-invariant spectrum (B(k) 2 ) = constant, the spikes 
are much less pronounced. This is expected , as the 
magnetic power k 2 (B(k) 2 } in that case is concentrated 
on small spatial scales, whereas for the Kolmogorov spec- 
trum, the power is concentrated on large scales. 



60. < E/EeV < 9000. 
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FIG. 4. The average over 4 magnetic field realizations of 
the distribution of events above 60 EeV in Supergalactic lat- 
itude for the scenario with B rms = 0.5 /iG shown in Figs. |l| 
andpnthick histogram), and for B rulB = 0.05 fi G, as shown in 
Figs^j] and |i] (thin histogram), assuming 1.6° angular resolu- 
tion. The dash-dotted curve represents a completely isotropic 
distribution. 

In Fig. [|, we show the histogram of events plotted vs. 
Supergalactic latitude; this figure clearly illustrates the 
relative anisotropy of a scenario with B rms — 0.05 /uG, 
with a distribution of sources strongly centered around 
Virgo, and the relative isotropy of a scenario with B lms = 
0.5 fiG. A similar plot for the observed AGASA events is 
given in Refs. M. 

Observational tests. A turbulent magnetic field 
of strength B lms ~ 0.5 fiG in the Local Super- 
cluster is not unrealistic [17y; it roughly corre- 



sponds to the expected equipartition value B cq ~ 

0.2 /jG (T/3 x 10 6 K) 1/2 (p b /0.3p c ) 1/2 , where T, p b , and 
p c respectively denote the temperature, baryon density 
and critical density jjj . The above scenario thus reopens 
a window for conventional astrophysical sources as ori- 
gin of UHECRs, provided that such sources can indeed 
accelerate protons up to E > 1000 EeV 0. In the follow- 
ing, we give several predictions and tests of the present 
scenario. 

The most direct test would be to detect and measure 
a permeating magnetic field in the Local Supercluster, 
or at least obtain a trace of its presence via the obser- 
vation of synchrotron emission. In this respect, we note 
the positive detection of a bridge of synchrotron emission 
in the plane of the Coma/Abell 1367 supercluster p8[ , 
which lends further support to our scenario. Unfortu- 
nately, synchrotron data does not yield the strength of 
the magnetic field, and measurements of Faraday rota- 
tion of polarized light from sources located in the Local 
Supercluster are needed to this end. 

Our model predicts large-scale anisotropy, whose am- 
plitude is expected to increase with energy above the 
pion production threshold, because of decreasing mag- 
netic deflection, and because sources located beyond the 
Supergalactic plane, that we neglected here, and that 
would tend to make the distribution more isotropic, do 
not contribute for E > 100 EeV |§. We quantify this 
anisotropy in terms of dipolc and quadrupole moments, 
assuming that there is full sky coverage, as planned for 
the Pierre Auger Observatory pQ] . For instance, the de- 
tection of 64 events above 60 EeV, corresponding to one 
year of one site of the Pierre Auger Observatory, would 
allow to detect anisotropy with a false alarm risk of 5% 
(resp. 1%) in 77. - 99.9% (resp. 63. - 99.7%) of cases 
if B rms = 0.05 /iG, and in 51. - 93% (resp. 26. - 83.%) 
of cases if B rms = 0.5 fj,G, where the ranges given cor- 
respond to variation of the source distribution, coher- 
ence length, Tig, and observer position. These statistics 
are limited by the finite number or events, not by cos- 
mic variance. We find that the anisotropy strongly de- 
creases with increasing radius of the source distribution, 
and with increasing magnetic field strength. Finally, we 
emphasize that a correlation of UHECR statistics with 
Supergalactic latitude is not a proper measure of the pre- 
dicted anisotropy when B Yms ~ 0.5 /iG, because this cor- 
relation starts to differ from the isotropic prediction only 
for E > 200 EeV (see also Fig. |). 

Finally, the observation of a repeated number of clus- 
ters of UHECRs with arrival directions compatible within 
the angular resolution of future instruments will also al- 
low to constrain severely the present model. A strong 
magnetic field leads to an almost isotropic arrival direc- 
tion distribution on large angular scales. At the same 
time, magnetic lensing, illustrated in Ref. j^] in the case 
of a point source, leads to significant small scale cluster- 
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ing. Overall, for any source distribution the probability 
for clustering of UHECRs in the presence of a magnetic 
field is higher than for an isotropic arrival direction dis- 
tribution, and (of course) smaller than in the absence of 
magnetic fields with point-like sources. 

TABLE I. Probabilities to detect 5 doublets above 40 EeV 
in the AGASA data set (47 events), and > l(resp.3) multi- 
plets of > 5(resp.3) events in a data set of 64 events above 
60 EeV, as expected for 1 year of the Southern Pierre Auger 
site (PAO), for B TU1S as indicated. 



-Brms (/iG) 


AGASA prob.(%) 


PAO prob.(%) 


(isotropic) 


0.29 


< 0.01 (< 0.01) 


0.05 


8. - 20. 


0.3-6. (0.5-4.) 


0.5 


8. - 16. 


5. - 8. (~ 11.) 



We quantify this in Table || by comparing cluster- 
ing probabilities for our scenario and for a completely 
isotropic distribution, in two cases: in the first, we give 
probabilities corresponding to the numbers of UHECR 
events and clusters observed by the AGASA experi- 
ment ||. To this end we simulated the finite angular 
resolution ~ 1.6° and sky coverage of the AGASA ex- 
periment, and calculated the probability of repeated oc- 
curences (multiplets). We call a multiplet a cluster of 
events such that all events fall within 2.5° of the first 
event detected. In the second case, we give an example 
for numbers of multiplets that would strongly discrim- 
inate our scenario against uniform source distributions 
with much weaker magnetic fields, for one year of the 
Southern Pierre Auger site, assuming 1° angular resolu- 
tion. The range of values given corresponds to variation 
of the source density profile, n^, coherence length, and 
position of the observer. 

These numbers indicate that the preference for strong 
fields increases with exposure. We also found that clus- 
tering increases with decreasing ns, with increasing ra- 
dius of the source density profile, and with the coherence 
length of the field. Within our model, a rather coherent 
field, with power on the largest scales, such as a Kol- 
mogorov spectrum is thus favored. Note that for a less 
diffuse source distribution with less abundant but more 
powerful sources, these probabilities would be larger. 
Our scenario, however, predicts the absence of correla- 
tions between the UHECR clusters and powerful sources 
associated with the large scale structure. This can be 
used to discriminate it also against models with highly 
structured source distributions and negligible magnetic 
fields plj| |. A somewhat larger number of clusters and/or 
higher multiplicity than in Tab. | would rule out our sce- 
nario, e.g., for a cluster of nine events out of 64 showers 
above 60 EeV, and B lms = 0.5 /iG, the confidence level is 
< 0.1%. 
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